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Abstract. We solve the following geometric problem, which arises in several three- 
dimensional applications in computational geometry: For which arrangements of two 
lines and two spheres in R 3 are there infinitely many lines simultaneously transversal to 
the two lines and tangent to the two spheres? 
■ We also treat a generalization of this problem to projective quadrics: Replacing the 

, spheres in R 3 by quadrics in projective space P 3 , and fixing the lines and one general 

quadric, we give the following complete geometric description of the set of (second) 
quadrics for which the 2 lines and 2 quadrics have infinitely many transversals and 
tangents: In the nine-dimensional projective space P 9 of quadrics, this is a curve of 
degree 24 consisting of 12 plane conies, a remarkably reducible variety. 



43 



X 



P 3 



> 

^J- . Introduction 

In ||14j| , one of us (Theobald) considered arrangements of k lines and A—k spheres in 
having infinitely many lines simultaneously transversal to the k lines and tangent to the 
CNl ■ 4— k spheres. Since for generic configurations of k lines and A—k spheres there are only 

finitely many common transversals/tangents, the goal was to characterize the non-generic 
configurations where the discrete and combinatorial nature of the problem is lost. One 
case left open was that of two lines and two spheres. We solve that here. 

A second purpose is to develop and present a variety of techniques from computational 
algebraic geometry for tackling problems of this kind. Since not all our readers are fa- 
miliar with these techniques, we explain and document these techniques, with the goal of 
increasing their applicability. For that reason, we first deal with the more general problem 
where we replace the spheres in R 3 by general quadratic surfaces (hereafter quadrics) in 
complex projective 3-space P 3 . In order to study the geometry of this problem, we fix two 
lines and a quadric in general position, and describe the set of (second) quadrics for which 
there are infinitely many common transversals/tangents in terms of an algebraic curve. 
It turns out that this set is an algebraic curve of degree 24 in the space P 9 of quadrics. 
Factoring the ideal of this curve shows that it is remarkably reducible: 

Theorem 1. Fix two skew lines t\ and £2 and a general quadric Q in P 3 . The closure of 
the set of quadrics Q' for which there are infinitely many lines simultaneously transversal 
to i\ and £2 and tangent to both Q and to Q' is a curve of degree 24 in the P 9 of quadrics. 
This curve consists of 12 plane conies. 
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We prove this theorem by investigating the ideal defining the algebraic curve describing 
the set of (second) quadrics. Based on this, we prove the theorem with the aid of a 
computer calculation in the computer algebra system Singular [|J. As explained in 
Section |3|, the success of that computation depends crucially on the preceeding analysis of 
the curve. Quite interestingly, there are real lines t\ and £2 and real quadrics Q such that 
all 12 components of the curve of second quadrics are real. In general, given real lines £\, 
£2, and a real quadric Q, not all of the 12 components are defined over the real numbers. 

While the beautiful and sophisticated geometry of our fundamental problem on lines 
and quadrics could be sufficient motivation to study this geometric problem, the original 
motivation came from algorithmic problems in computational geometry. As explained 
in fl4 |, problems of this type occur in applications where one is looking for a line or ray 



interacting (in the sense of "intersecting" or in the sense of "not intersecting" ) with a given 
set of three-dimensional bodies, if the class of admissible bodies consists of polytopes and 
spheres (respectively quadrics). Concrete application classes of this type include visibility 
computations with moving viewpoints |15| , controlling a laser beam in manufacturing |TTJ , 



or the design of envelope data structures supporting ray shooting queries (i.e., seeking the 
first sphere, if any, met by a query ray) JO. With regard to related treatments of the 



resulting algebraic-geometric core problems, we refer to || [T(], In these papers, the 
question of arrangements of four (unit) spheres in M 3 leading to an infinite number of 
common tangent lines is discussed from various viewpoints. 

The present paper is structured as follows: In Section [I], we review the well-known 
Plucker coordinates from line geometry. In Section |2|, we characterize the set of lines 
transversal to two skew lines and tangent to a quadric in terms of algebraic curves; we 
study and classify these so-called (2, 2)-curves. Then, in Section |3], we study the set of 
quadrics which (for prescribed lines £\ and £ 2 ) lead to most (2, 2)-curves. This includes 
computer-algebraic calculations, based on which we establish the proof of Theorem [TJ. The 
appendix to the paper contains annotated computer code used in the proof. In Section £|, 
we give some detailed examples illustrating the geometry described by Theorem |], and 
complete its proof. Finally, in Section |5|, we solve the original question of spheres and 
give the complete characterization of configurations of two lines and two quadrics having 
infinitely many lines transversal to the lines and tangent to the quadrics. For a precise 



statement of that characterization see Theorems 14 and 18 



1. Plucker Coordinates 

We review the well-known Plucker coordinates of lines in three-dimensional (complex) 
projective space P 3 . For a general reference, see [[7], O. Let x = (x , X\, x 2 , x 3 ) T and 
V — iyoi Vii 2/2, yz) T £ P 3 be two points spanning a line £. Then £ can be represented (not 
uniquely) by the 4 x 2-matrix L whose two columns are x and y. The Plucker vector 
P — (Poi,Po2,Po3,Pi2,Pi3,P23) T G P 5 of £ is defined by the determinants of the 2x2- 
submatrices of L, that is, pij := XiUj — Xj]ji. The set Gj.^ of all lines in P 3 is called the 
Grassmannian of lines in P 3 . The set of vectors in P 5 satisfying the Plucker relation 



(1.1) 



P01P23 - P02P13 + P03P12 = 
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is in 1-1-correspondence with Gi^. See, for example Theorem 11 in § 8.6 of 0]. 

A line £ intersects a line £' in P 3 if and only if their Pliicker vectors p and p' satisfy 



;i-2) 



P01P23 - P02K3 + PosPu + P12P03 - P13P02 + PaaPoi = °- 



Geometrically, this means that the set of lines intersecting a given line is described by a 
hyperplane section of the Pliicker quadric ( |1.1| ) in P 5 . 

In Pliicker coordinates we also obtain a nice characterization (given in P3"| ) of the 
lines tangent to a given quadric in P 3 . (See JR] for an alternative deduction of that 
characterization). We identify a quadric x T Qx = in P 3 with its symmetric 4x4- 
representation matrix Q. Thus the sphere with center (ci,c 2 ,Cs) T G M 3 and radius r and 



described in P 3 by 
matrix 



£1 - cix ) 2 + (x 2 - c 2 x ) 2 + (x 3 - c 3 x y 



t 2 Xq, is identified with the 



{ci + 4 + 4 



\ 
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The quadric is smooth if its representation matrix has rank 4. To characterize the tangent 
lines, we use the second exterior power of matrices 



dm 
2 



(see [12 , p. 145l)||13|). Here C axb is the set of a x b matrices with complex entries. The 
row and column indices of the resulting matrix are subsets of cardinality 2 of {1, ... , m} 
and {1, . . . , n}, respectively. For / = {11,12} with 1 < i\ < i 2 < m and J = {ji,j 2 } with 
1 < h < h < n, 



A. .A. . 
^11,31^12,32 



A- -A- ■ 

^1,32^2,31 



Let i be a line in P 3 and L be a 4 x 2-matrix representing £. Interpreting the 6 x 1-matrix 
A 2 L as a vector in P 5 , we observe that A 2 L = p#, where pi is the Pliicker vector of £. 

Recall the following algebraic characterization of tangency: The restriction of the qua- 
dratic form to the line £ is singular, in that either it has a double root, or it vanishes 
identically. When the quadric is smooth, this implies that the line is tangent to the 
quadric in the usual geometric sense. 



Proposition 2 (Proposition 5.2 of [13|). A line £ C P 3 is tangent to a quadric Q if and 
only if the Pliicker vector pi of £ lies on the quadratic hypersurface in P 5 defined by A 2 Q, 
if and only if 



(1.3) 



pj (A 2 Q) Pi = 0. 
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For a sphere with radius r and center (ci,c 2 ,c 3 ) T G M. 3 the quadratic form pj \/\ 2 Q)pe. 

is 



(1.4) 
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T 


/ 4 


+ 4- 


r 2 -CiC 2 


-C1C3 


c 2 


c 3 





\ 


/ Poi \ 


P02 






-C\C 2 


c? + 4 - 


r 2 -c 2 c 3 


-ci 





c 3 




P02 


P03 






-C1C3 


-c 2 c 3 


c i + c 2 — r2 





-c\ 


-c 2 




P03 


Pl2 






c 2 


-ci 





1 










Pl2 


P13 






c 3 





-ci 





1 







Pia 


V P23 / 




\ 





c 3 


-c 2 








1 


/ 


\ P23 1 



2. Lines in P 3 meeting 2 lines and tangent to a quadric 

We work here over the ground field C. First suppose that £\ and £2 are lines in P 3 that 
meet at a point p and thus span a plane II. Then the common transversals to £\ and £ 2 
either contain p or they lie in the plane II. This reduces any problem involving common 
transversals to £\ and £2 to a planar problem in P 2 (or M 2 ), and so we shall always assume 
that i\ and £2 are skew. Such lines have the form 

£1 = {wa + xb : \w, x] G P 1 } , 

(2.1) 1 

£2 = {yc + zd : [^JeP 1 } 

where the points a, b,c,d G P 3 are affinely independent. We describe the set of lines 
meeting l\ and £ 2 that are also tangent to a smooth quadric Q. We will refer to this set 
as the envelope of common transversals and tangents, or (when £1 and £2 are understood) 
simply as the envelope of Q. 

The parameterization of Q2.1Q allows us to identify each of £\ and £2 with P 1 ; the point 
wa + xb G £\ is identified with the parameter value [w,x] G P 1 , and the same for £ 2 . We 
will use these identifications throughout this section. In this way, any line meeting l\ and 
£2 can be identified with the pair ([w, x], [y, z\) G P 1 x P 1 corresponding to its intersections 
with l\ and £ 2 - By ( |1.2|) , the Pliicker coordinates pe = pe(w, x,y, z) of the transversal £ 
passing through the points wa + xb and yc + zd are separately homogeneous of degree 1 in 
each set of variables {w, x} and {y,z}, called bihomogeneous of bidegree (1,1) (see, e.g., 
I §H..Y). 

By Proposition |2|, the envelope of common transversals to £\ and £2 that are also tangent 
to Q is given by the common transversals £ of £\ and l 2 whose Pliicker coordinates pi 
additionally satisfy pAf\ 2 Q)pe. = 0. This yields a homogeneous equation 

(2.2) F(w,x,y,z) := p e (w, x,y, z) J (A 2 Q)p e (w 1 x,y, z) = 

of degree four in the variables w, x, y, z. More precisely, F has the form 

2 

(2.3) F(w,x,y,z) = £ c^x^Y^ 

i,j=0 

with coefficients c^-, that is F is bihomogeneous with bidegree (2, 2). The zero set of a (non- 
zero) bihomogeneous polynomial defines an algebraic curve in P 1 x P 1 (see the treatment of 
projective elimination theory in 0, §8.5]). In correspondence with its bidegree, the curve 
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defined by F is called a (2, 2)-curve. The nine coefficients of this polynomial identify the 
set of (2, 2)-curves with P 8 . 

It is well-known that the Cartesian product P 1 x P 1 is isomorphic to a smooth quadric 
surface in P 3 Proposition 10 in § 8.6]. Thus the set of lines meeting t\ and £2 and 
tangent to the quadric Q is described as the intersection of two quadrics in a projective 3- 
space. When it is smooth, this set is a genus 1 curve || Exer. 1.7.2(d) and Exer. 11.8.4(g)]. 
This set of lines cannot be parameterized by polynomials — only genus curves (also called 
rational curves) admit such parameterizations (see, e.g., [§, Corollary 2 on p. 268]). This 
observation is the starting point for our study of common transversals and tangents. 

Let C be a (2, 2)-curve in P 1 xP 1 defined by a bihomogeneous polynomial F of bidegree 2. 
The components of C correspond to the irreducible factors of F, which are bihomogeneous 
of bidegree at most (2, 2). Thus any factors of F must have bidegree one of (2, 2), (2, 1), 
(1, 1), (1, 0), or (0, 1). (Since we are working over C, a homogeneous quadratic of bidegree 
(2,0) factors into two linear factors of bidegree (1,0).) Recall (for example, 0) a point 
([wo,^o]; [2/0 ; 20]) € C C P 1 x P 1 is singular if the gradient VF vanishes at that point, 
VF([w ,x ], [2/0 j 20]) = 0. The curve C is smooth if it does not contain a singular point; 
otherwise C is singular. We classify (2, 2)-curves, up to change of coordinates on l\ x £ 2 , 
and interchange of £\ and £2- Note that an (a, 6)-curve and a (c, d)-curve meet if either 
ad 7^ or be 7^ 0, and the intersection points are singular on the union of the two curves. 

Lemma 3. Let C be a (2,2)-curve on P 1 x P 1 . Then, up to interchanging the factors of 
P 1 x P 1 , C 2S either 

1. smooth and irreducible, 

2. singular and irreducible, 

3. the union of a (1,0) -curve and an irreducible (1,2) -curve, 

4. the union of two distinct irreducible (1, l)-curves, 

5. a single irreducible (1,1) -curve, of multiplicity two, 

6. the union of one irreducible (1,1) -curve, one (l,0)-curve, and one (0,1) -curve, 

7. the union of two distinct (l,0)-curves, and two distinct (0, l)-curves, 

8. the union of two distinct (1,0) -curves, and one (0,1) -curve of multiplicity two, 

9. the union of one (l,0)-curve, and one (0,l)-curve, both of multiplicity two. 
In particular, when C is smooth it is also irreducible. 

When the polynomial F has repeated factors, we are in cases (5), (8), or (9). We study 
the form F when the quadric is reducible, that is either when Q has rank 1, so that it 
defines a double plane, or when Q has rank 2 so that it defines the union of two planes. 

Lemma 4. Suppose Q is a reducible quadric. 

(1) If Q has rank 1, then A 2 Q = 0, and so the form F in Q2.2| ) is identically zero. 

(2) Suppose Q has rank 2, so that it defines the union of two planes meeting in a line 
i. If i is one of £1 or £2, then the form F in ( j2.2|) is identically zero. Otherwise 
the form F is the square of a (1, l)-form, and hence we are in cases (5) or (9) of 
Lemma [|. 

Proof. The first statement is immediate. For the second, let £' be a line in P 3 with 
Plucker coordinates p^. From the algebraic characterization of tangency of Proposition |2|, 
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pJi(/\ 2 Q)pe> = implies that the restriction of the quadratic form to £' either has a zero 
of multiplicity two, or it vanishes identically. In either case, this implies that £' meets the 
line £ common to the two planes. Conversely, if £' meets the line £, then pJi(A 2 Q)pi' = 0. 

Thus if £ equals one of £\ or £2, then pJi(A 2 Q)pe/ = for every common transversal £' 
to £\ and £2, and so the form F is identically zero. Suppose that £ is distinct from both 
£\ and £2- We observed earlier that the set of lines transversal to £\ and £2 that also meet 
£ is defined by a (1, l)-form G. Since the (2, 2)-form F defines the same set as does the 
(1, l)-form G, we must have that F = G 2 , up to a constant factor. □ 

As above, let C be defined by the polynomial F. For a fixed point [w, x], the restriction 
of the polynomial F to [w,x] x P 1 is a homogeneous quadratic polynomial in y,z. A 
line passing through [w,x] £ £\ and the point of £2 corresponding to any zero of this 
restriction is tangent to Q. This construction gives all lines tangent to Q that contain the 
point [w, x\. We call the zeroes of this restriction the fiber over [w, x] of the projection of 

C to /;. 

We investigate these fibers. Consider the polynomial F as a polynomial in the variables 
y, z with coefficients polynomials in w,x. The resulting quadratic polynomial in y, z has 
discriminant 

(2.4) f^ca^V-A - i(j2^x 2 A (j^ ^ 2 '' 

\i=0 / \i=0 J \i=0 

Lemma 5. If this discriminant vanishes identically, then the polynomial F has a repeated 
factor. 

Proof. Let a, (3, 7 be the coefficients of y 2 ,yz, z 2 in the polynomial F, respectively. Then 
we have (3 2 = 4a7, as the discriminant vanishes. Since the ring of polynomials in w, x is 
a unique factorization domain, either a differs from 7 by a constant factor, or else both 
a and 7 are squares. If a and 7 differ by a constant factor, then so do a and (3. Writing 
(3 = Ida for some d £ C, we have 

F = ay 2 + Idayz + d 2 az 2 = a(y + dz) 2 . 

If we have a = 5 2 and 7 = a 2 for some linear polynomials 5 and cr, then 

F = 5 2 y±25ayz + a 2 z 2 = (Sy±az) 2 . 

□ 

When F does not have repeated factors, the discriminant does not vanish identically. 
Then the fiber of C over the point [w, x] of £\ consists of two distinct points exactly when 
the discriminant does not vanish at Since the discriminant has degree 4, there are 

at most four fibers of C consisting of a double point rather than two distinct points. We 
call the points [w,x] of £\ whose fibers consist of such double points ramification points 
of the projection from C to £\. 

This discussion shows how we may parameterize the curve C, at least locally. Suppose 
that we have a point [w,x] £ P 1 where the discriminant G2.4Q does not vanish. Then we 
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may solve for [y,z] in the polynomial F in terms of [u>,x]. The different branches of the 
square root function give local parameterizations of the curve C. 

2.1. A normal form for asymmetric smooth (2, 2)-curves. Recall that for any dis- 
tinct points ax, a,2,a3 G P 1 and any distinct points 61,62,^3 €= P 1 ? there exists a pro- 
jective linear transformation (given by a regular 2 x 2-matrix) which maps a« to bi, 
l<i<3||. 

Lemma 6. If the (2, 2) -curve is smooth then the projection of C to t\ has four different 
ramification points. 

Proof. Changing coordinates on l\ and £2 by a projective linear transformation if neces- 
sary, we may assume that this projection to t\ is ramified over [w,x] = [1,0], and the 
double root of the fiber is at [y, z] = [1, 0]. Restricting the polynomial F ( |2.3|) to the fiber 
over [w,x] = [1,0] gives the equation 

C22V 2 + oiiyz + c 20 z 2 = 0. 

Since we assumed that this has a double root at [y, z] = [1, 0], we have C21 = C22 = 0. 

Suppose now that the projection from C to i\ is ramified at fewer than four points. We 
may assume that [w,x] = [1,0] is a double root of the discriminant (|2.4|) , which implies 
that the coefficients of w A and w 3 x in ( |2.4| ) vanish. The previously derived condition 
C21 = C22 = implies that the coefficient of w A vanishes and the coefficient of w 3 x becomes 
— 4c 2 oCi2- If c 2 o = 0, then every non-vanishing term of ( |2.3| ) depends on x; hence, x divides 
F, and so C is reducible, and hence not smooth. If c±2 = then the gradient VF vanishes 
at the point ([1,0], [1,0]), and so C is not smooth. □ 

Suppose that C is a smooth (2, 2)-curve. Then its projection to t\ is ramified at four 
different points. We further assume that the double points in the ramified fibers project 
to at least 3 distinct points in £ 2 - We call such a smooth (2, 2)-curve asymmetric. The 
choice of this terminology will become clear in Section |4j. We will give a normal form for 
such asymmetric smooth curves. 

Hence, we may assume that three of the ramification points are [w,x] = [0,1], [1,0], 
and [1, 1], and the double points in these ramification fibers occur at [y, z] = [0, 1], [1, 0], 
and [1,1], respectively. As in the proof of Lemma §, the double point at [y,z] = [1,0] 
in the fiber over [w,x] = [1,0] implies that C21 = C22 = 0. Similarly, the double point 
at [y, z] = [0, 1] in the fiber over [w, x] = [0, 1] implies that c 00 = c i = 0. Thus the 
polynomial F (|2.3|) becomes 

C2qW 2 z 2 + CiqWxz 2 + Cnwxyz + C\2wxy 2 + CQ 2 x 2 y 2 
Restricting F to the fiber of [w,x] = [1, 1] gives 

c w z 2 + c 20 z 2 + c u yz + c 02 y 2 + c 12 y 2 . 
Since this has a double root at [y, z] = [1,1], we must have 

1 

— -Cn — C10 + C20 — C02 + C\2 ■ 
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Dehomogenizing (setting en = 
obtain the following theorem. 



-2) and letting c 2 q := s and C02 := t for some s, t G C, we 



Theorem 7. After projective linear transformations in £1 and £ 2 , an asymmetric smooth 
(2, 2)-curve is the zero set of a polynomial 



sw 2 z 2 + (1- 



■s )wxz 



— 2wxyz + (l—t)wxy" 



+ tx 2 y 2 . 



of the polynomial (|2.5| ) is 



(2.5) 

for some (s, t) G C 2 satisfying 

(2.6) st(s-l)(f-l)(s-*) ^ 0. 

We complete the proof of Theorem |7|. The discriminant 

4wx(w—x) (s(t—l)w — t(s-l)x) , 

which has roots at [w,x] = [0, 1], [1, 0], [1, 1], and a = [t(s — 1), s(t — 1)]. Since we assumed 
that these are distinct, the fourth point a must differ from the first three, which implies 
that (s, t) satisfies (2.6). The double point in the fiber over a occurs at [y, z] = [s — 1, t— 1]. 
This equals a double point in another ramification fiber only for values of the parameters 
not allowed by (127 



Remark 8. These calculations show that smooth (2, 2)-curves exhibit the following di- 
chotomy. Either the double points in the ramification fibers project to four distinct points 
in £ 2 or to two distinct points. They must project to at least two points, as there are at 
most two points in each fiber of the projection to £2- We showed that if they project to 
at least three, then they project to four. 

We compute the parameters s and t from the intrinsic geometry of the curve C. Recall 
the following definition of the cross ratio (see, for example [0, §1.1.4]). 

G P 1 with ai = [ai,Pi\, the cross ratio of 



Definition 9. For four points a%, . . . ,04 
ax, . . . ,04 is the point of P 1 defined by 



det 



det 



a.2 0:4 
fa A 



det 



Pi 



«3 
ft 



det 



a 2 

02 



«3 
03 



If the points are of the form 



this simplifies to 

04-0! 4 ~02 



_03 ~ 01 ' 03 ~02_ 

The cross ratio of four points 01,02,03,04 G P 1 remains invariant under any projective 
linear transformation. 

The projection of C to £\ is ramified over the points [w,x] = [0, 1], [1, 0], [1, 1] and 
a = [t(s — 1), s(t — 1)]. The cross ratio of these four (ordered) ramification points is 
[t(s— 1), s(t— 1)]. Similarly, the cross ratio of the four (ordered) double points in the 
ramification fibers is [s— 1, t—1]. 
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This computation of cross ratios allows us to compute the normal form of an asymmetric 
smooth (2, 2)-curve. Namely, let 01,02,03, and 04 be the four ramification points of the 
projection of C to t\ and 61, b 2 , 63, and 64 be the images in £ 2 of the corresponding double 
points. Let 71 be the cross ratio of the four points a 1; a 2 , a 3 , and a 4 (this is well-defined, 
as cross ratios are invariant under projective linear transformation). Similarly, let 72 be 
the cross ratio of the points bi,b 2 ,b 3 , and 64. For four distinct points, the cross ratio is 
an element of C \ {0, 1}, so we express 71,72 as complex numbers. The invariance of the 
cross ratios yields the conditions on s and t 

s(t-l) , t-1 

- — = 71 and — - = 73. 

t[s— 1) s—1 

Again, since 71, 72 G C \ {0, 1}, these two equations have the unique solution 

_ 7i(72-l) , _ 72-1 

b — ~ — dllLl i — . 

72(71 - 1) 71 - 1 

3. Proof of Theorem 1 

We characterize the quadrics Q which generate the same envelope of tangents as a given 
quadric. A symmetric 4x4 matrix has 10 independent entries which identifies the space 
of quadrics with P 9 . Central to our analysis is a map <f defined for almost all quadrics 
Q. For a quadric Q (considered as a point in P 9 ) whose associated (2, 2)-form ( |2.2p is not 
identically zero, we let <f(Q) be this (2, 2)-form, considered as a point in P 8 . With this 
definition, we see that the Theorem |l] is concerned with the fiber <y9 -1 (C), where C is the 
(2, 2)-curve associated to a general quadric Q. Since the domain of <f is 9-dimensional 
while its range is 8-dimensional, we expect each fiber to be 1-dimensional. 

We will show that every smooth (2, 2) curve arises as f{Q) for some quadric Q. It is 
these quadrics that we meant by general quadrics in the statement of Theorem p]. This 
implies that Theorem [l] is a consequence of the following theorem. 



Theorem 10. Let C G P 8 be a smooth (2, 2) -curve. Then the closure y>~ 1 (C) in P 9 of 
the fiber of if is a curve of degree 24 that is the union of 12 plane conies. 

We prove Theorem [l^ by computing the ideal J of the fiber if {C). Then we factor 
J into several ideals, which corresponds to decomposing the curve of degree 24 into the 
union of several curves. Finally, we analyze the output of these computations by hand to 
prove the desired result. 

Our initial formulation of the problem gives an ideal / that not only defines the fiber 
of if, but also the subset of P 9 where ip is not defined. We identify and remove this subset 
from / in several costly auxiliary computations that are performed in the computer algebra 
system Singular M. It is only after removing the excess components that we obtain the 
ideal J of the fiber <^ -1 (C). 

Since we want to analyze this decomposition for every smooth (2, 2)-curve, we must 
treat the representation of C as symbolic parameters. This leads to additional difficulties, 
which we circumvent. It is quite remarkable that the computer-algebraic calculation 
succeeds and that it is still possible to analyze its result. 
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(3.3) Q 



In the following, we assume that l\ is the x-axis. Furthermore, we may apply a projec- 
tive linear transformation and assume without loss of generality that £2 is the yz-Yme at 
infinity. Thus we have 

£ 1 = {(w, x, 0, 0) T G P 3 : [w, x] G P 1 } , 
£2 = {(0,0,y,z) T GP 3 : MGP 1 }. 
Hence, in Pliicker coordinates, the lines intersecting l\ and £2 are given by 

(3.1) {(0, wy, wz, xy, xz, 0) T G P 5 : [w, x], [y, z] G P 1 } . 

By Proposition ^ the envelope of common transversals to £\ and £2 that are also tangent 
to Q is given by those lines in fl3.1| ) which additionally satisfy 

(3.2) (0, wy, wz, xy, xz, 0) (A 2 Q) (0,wy,wz,xy,xz,0) T = 0. 

A quadric Q in P 3 is given by the quadratic form associated to a symmetric 4 x 4-matrix 

/ a b c d\ 
b e f g 
c f h k 
\d g k I J 

In a straightforward approach the ideal / of quadrics giving a general (2, 2)-curve C is 
obtained by first expanding the left hand side of (|3.2|) into 

(el-g 2 )x 2 z 2 + 2{bl-dg)wxz 2 + (al-d?)w 2 z 2 

(3.4) + 2(ek-gf)x 2 yz + 2(2bk-cg-df)wxyz + 2(ak-dc)w 2 yz 

+ (eh—f 2 )x 2 y 2 + 2{bh—cf)wxy 2 + (ah—c 2 )w 2 y 2 . 

We equate this (2, 2) -form with the general (2, 2)-form (|2.3|) , as points in P 8 . This is 
accomplished by requiring that they are proportional, or rather that the 2x9 matrix of 
their coefficients 

coo C10 c 20 c i c u ■■■ C 2 2 

el-g 2 2{bl-dg) al - d 2 2(ek-gf) 2(2bk - eg - df) ... ah - c 2 

has rank 1. Thus the ideal / is generated by the (2) minors of this coefficient matrix. 

With this formulation, the ideal / will define the fiber (^ _1 (C) as well as additional, 
excess components that we wish to exclude. For example, the variety in P 9 defined by 
the vanishing of the entries in the second row of this matrix will lie in the variety /, 
but these points are not those that we seek. Geometrically, these excess components 
are precisely where the map tp is not defined. By Lemma [|, we can identify three of 
these excess components, those points of P 9 corresponding to rank 1 quadrics, and those 
corresponding to rank 2 quadrics consisting of the union of two planes meeting in either £\ 
or in £2. The rank one quadrics have ideal E\ generated by the entries of the matrix A 2 Q, 
the rank 2 quadrics whose planes meet in £1 have ideal E 2 generated by a, b, c, d, e, f, g, 
and those whose plane meets in £ 2 have ideal E 3 generated by c, d, f, g, h, k, I. 

We remove these excess components from our ideal I to obtain an ideal J whose set of 
zeroes contain the fiber y9 _1 (C). After factoring J into its irreducible components, we will 
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observe that (f does not vanish identically on any component of J, completing the proof 
that J is the ideal of (^ _1 (C), and also the proof of Theorem [H]. 

Since coo, cio, . . . , C22 have to be treated as parameters, the computation should be 
carried out over the function field Q(c 00 , Cxo, • • ■ , C22). That computation is infeasible. 
Even the initial computation of a Grobner basis for the ideal I (a necessary prerequisite) 
did not terminate in two days. In contrast, the computation we finally describe termi- 
nates in 7 minutes on the same computer. This is because the original computation in 
Q( c oo, cio, • • • , C22)[a, b, . . . , /] involved too many parameters. 



We instead use the 2-parameter normal form (|2.5|) for asymmetric smooth (2, 2)-curves. 



This will prove Theorem 1C in the case when C is an asymmetric smooth (2, 2)-curve. We 
treat the remaining cases of symmetric smooth (2, 2)-curves in Section [|. As described in 
Section [Oj , by changing the coordinates on i\ and £2, every asymmetric smooth (2,2)- 
curve can be transformed into one defined by a polynomial in the family (|2.5|) . Equating 
the (2, 2)-form ( |3.4|) with the form (|2.5|) gives the ideal I generated by the following 
polynomials: 

(3.5) el — g 2 , ek — gf, ak — dc, ah — c 2 , 

and the ten 2x2 minors of the coefficient matrix: 

(o fi \ M - ( s 1-5 ~ 2 1 ~ 1 t 

■ \al-d 2 2{bl - dg) 2{2bk -eg- df) 2(bh - cf) eh - f 2 

This ideal I defines the same three excess components as before, and we must remove 
them to obtain the desired ideal J. Although the ideal I should be treated in the ring 
S := Q(s, t)[a, b, c, d, e, f, g, h, k, I], the necessary calculations are infeasible even in this 
ring, and we instead work in subring R := Q[a, b, c, d, e, f, g, h, k, l][s, £]. In the ring R, 
the ideal I is homogeneous in the set of variables a,b, . . . , / , thus defining a subvariety 
of P 9 x C 2 . The ideals Ei, E 2 , and E 3 describing the excess components satisfy Ej D J, 
1<J<3. 

A Singular computation shows that J is a five-dimensional subvariety of P 9 x C 2 (see 
the Appendix for details). Moreover, the dimensions of the three excess components are 
5, 4, and 4, respectively. In fact, it is quite easy to see that dim E 2 = dim E 3 = 4 as 
both ideals are defined by 7 independent linear equations. 

We are faced with a geometric situation of the following form. We have an ideal I whose 
variety contains an excess component defined by an ideal E and we want to compute the 
ideal of the difference 

V(/) - V(E) , 

here, V(K) is the variety of an ideal K. Computational algebraic geometry gives us an 
effective method to accomplish this, namely saturation. The elementary notion is that of 
the ideal quotient (J : E), which is defined by 

(I:E) := {feR\fgeI for all g G E} . 
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Then the saturation of I with respect to E is 

oo 

(7 : E°°) : = |J(J: E n ) . 

n=l 

The least number n such that (/ : E°°) = (I : E n ) is called the saturation exponent. 

Proposition 11 (0, §4.4] or |3|, §15.10] or the reference manual for Singular). Over an 
algebraically closed field, 

V(I:E°°) = V(I) - V(E) . 

A Singular computation shows that the saturation exponent of the first excess ideal 
Ei in I is 1, and so the ideal quotient suffices to remove the excess component V(-E'i) from 
V(i). Set I' := (J : Ei), an ideal of dimension 4. The excess ideals E 2 and E 3 each have 
saturation exponent 4 in Ii, and so we saturate /' with respect to each to obtain an ideal 
J := ( (/' : E%°) : E%°), which has dimension 3 in P 9 x C 2 . 

To study the components of V(J), we first apply the factorization Grobner basis al- 
gorithm to J, as implemented in the Singular command f acstd (see [|5| or the refer- 
ence manual of Singular). This algorithm takes two arguments, an ideal I and a list 
L = fx, . . . , f n of polynomials. It proceeds as in the usual Buchberger algorithm to com- 
pute a Grobner basis for J, except that whenever it computes a Grobner basis element G 
that it can factor, it splits the calculation into subcalculations, one for each factor of G 
that is not in the list L, adding that factor to the Grobner basis for the corresponding 
subcalculation. The output of f acstd is a list Jj, 1%, . . . , I m of ideals with the property 
that 



\JV(Ij) -V(ff-f n ) = V(J) - V(/i ••■/„) . 

Thus, the zero set of I coincides with the union of zero sets of the factors Ij, in the region 
where none of the polynomials in the list L vanish. In terms of saturation, this is 

(3.7) rad(/ 1 • • • I m : (fj 2 ■ ■ ■ f n )°°) = rad(J : (fj 2 ■ ■ ■ f n )°°) 

where rad(if) denotes the radical of an ideal K. Some of the ideals Ij may be spurious 
in that V(Ij) is already contained in the union of the other V(i»). 

We run f acstd on the ideal J with the list of polynomials s, t, s—1, t—1, and s—t, and 
obtain seven components J , Ji, . . . , Jq. The components Ji, . . . , J 6 each have dimension 
3, while the component Jq has dimension 2. Since V(Jq) is contained in the union of the 
V(Ji), . . . , V(J 6 ), it is spurious and so we disregard it. 

We now, finally, change from the base ring R to the base ring S, and compute with 
the parameters s, t. There, J defines an ideal of dimension 1 and degree 24 in the 9- 
dimensional projective space over the field Q(s,t). As we remarked before, we have that 
V( J) D V9 _1 (C). The factorization of J into J\, . . . , Jq remains valid over S. The reason 
we did not compute the factorization over S is that f acstd and the saturations were 
infeasible over S, and the standard arguments from computational algebraic geometry we 
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have given show that it suffices to compute without parameters, as long as care is taken 
when interpreting the output. 

Each of the factors Jj has dimension 1 and degree 4. Moreover, each ideal contains 
a homogeneous quadratic polynomial in the variables k, I which must factor over some 
field extension of Q(s,i). In fact, these six quadratic polynomials all factor over the field 
Q(\/s, \fi). For example, two of the Jj contain the polynomial (s — l)k 2 — 2kl — I 2 , which 
is the product 

((Vi+l)* + ((Vi-l)*; - • 
For each ideal Jj, the factorization of the quadratic polynomial induces a factorization of 
Ji into two ideals Jn and J^. Inspecting a Grobner basis for each ideal shows that each 
defines a plane conic in P 9 . Thus, over the field Q(y/s, y/i), J defines 12 plane conies. 
Theorem [TTJ is a consequence of the following two observations. 

(1) The factorization of J gives 12 distinct components for all values of the parameters 
s, t satisfying ( |2.6| ). 

(2) The map <p does not vanish identically on any of the components V(Jij) for values 
of the parameters s,t satisfying (|2.6|). 

By (1), no component of J is empty for any s,t satisfying ( j2.6|) and thus, for every 
asymmetric (2,2)-curve C, there is a quadric Q with ip(Q) = C. Also by (1), J has 
exactly 12 components with each a plane conic, for any s,t satisfying fl2.6|) , and by (2), 

v( J) = VW)- 

4. Symmetric smooth (2,2)-curves 

We investigate smooth curves C whose double points in the ramified fibers over i\ 
have only two distinct projections to £2- Assume that the ramification is at the points 
[w, x] = [1, 1], [1, —1], [1, s], and at [1, — s], for some s G C\ {0, ±1} with the double points 
in the fibers at [y, z] = [1, 0] for the first two and at [0, 1] for the second two. Since the 
points [1, 1], [1, —1], [1, s], and [1, — s] have cross ratio 













[ ' (l + s) 2 J 



we see that all cross ratios in P 1 \{[1, 0], [0, 1], [1, 1]} are obtained for some s G C\{0, ±1}. 
Thus our choice of ramification results in no loss of generality. 

As in Section 0, these conditions give equations on the coefficients Cy of the general 
(2,2)-curve (p|): 

coo + cio + c 20 = 0, c i + c u + C21 = 0, coo - cio + c 20 = 0, 
coi - cu + c 2 i = 0, c 02 + ci 2 s + c 22 s 2 = 0, cqi + ens + C21S 2 = 0, 

c 02 — c 12 s + C22S 2 = 0, Cqi — C\\S + C2lS 2 = 0. 

These equations have the following consequences 

= C21 = Cqi = C\2 = C\\ = Cio = C02 + C22S 2 = Coo + c 20 • 

Hence after normalizing by setting C20 = 1, the (2,2)-form ( ^.3| ) becomes 

/ 2 2\ 2 1 I 2 2 2\ 2 

[X — W )y + C22{X — s w )z . 
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While the choice of ramification points [1, 1], [1, —1], [1, s], [1, — s] fixes the parameteriza- 
tion of £i, the double points in the fibers of [1, 0] and [0, 1] do not fix the parameterization 
of £2- Thus we are still free to scale the ^-coordinate. We normalize this equation setting 
C22 = ±1. We do not simply set C22 = 1 because that misses an important real form of 
the polynomial. This normalization gives 

(4.1) (x 2 -w 2 )y 2 ± (x 2 -s 2 w 2 )z 2 = (y 2 ±z 2 )x 2 -(y 2 ±s 2 z 2 )w 2 . 

This shows the equation to be symmetric under the involution [w,x] <-> [V^Tz, y\- This 
symmetry is the source of our terminology for the two classes of (2, 2)-curves. Also, if 
s ^ {±1,0}, then this is the equation of a smooth (2, 2)-curve. With the choice of sign 
(— ), which we call the curve C(s). 

Note that (|4. 1|) is real if s either is real or is purely imaginary (s £ M.y/—1 ). We 
complete the proof of Theorem [1] with the following result for symmetric (2, 2)-curves. 



Theorem 12. For each s £ C \ {±1, 0} ; the closure of the fiber <yj _1 (C(s)) consists of 12 
distinct plane conies. When s £ K or s £ M.y/—L and we use the real form of ( |4. 1| ) with 
the plus sign (+), then exactly 4 of these 12 components will be real. If we use the real 
form of ( |4.1| ) with the minus sign {—), then if s £ R, all 12 components will be real, but 
if s £ 1; then exactly 4 of these 12 components will be real. 

Proof. Our proof follows the proof of Theorem ^ almost exactly, but with significant 
simplifications and a case analysis. Unlike the proof described in Section we do not give 
annotated Singular code in an appendix, but rather supply such annotated Singular 
code on the web page^. 

The outline is as before, except that we work over the ring of parameters Q(s), and 
find no extraneous components when we factor the ideal into components. We formulate 
this as a system of equations, remove the same three excess components, and then factor 
the resulting ideal. We do this calculation four times, once for each choice of sign (±) 



in (|4.1|) , and for s £ M and s £ 1- Examining the output proves the result. □ 



We consider in some detail four cases of the geometry studied in Section 0, which 



correspond to the four real cases of Theorem 12. As in Section |2|, let £\ be the x-axis 
and I2 be the yz-\ine at infinity. Viewed in M 3 , lines transversal to £\ and I2 are the set 
of lines perpendicular to the x-axis. For a transversal line £, the coordinates [y, z] of the 
point £ fl £2 can be interpreted as the slope of £ in the two-dimensional plane orthogonal 
to the x-axis. 

Consider real quadrics given by an equation of the form 
(4.2) x 2 + (y-y ) 2 ±z 2 = 1. 

The quadrics with the plus (+) sign are spheres with center (0,y ,0) T and radius 1, 
and those with the minus (— ) sign are hyperboloids of one sheet. When \yo\ > 1 the 
quadric does not meet the x-axis. We look at four families of such quadrics: spheres and 
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hyperboloids that meet and do not meet the x-axis. We remark that quadrics which are 
tangent to the x-axis give singular (2, 2)-curves. 
First, consider the resulting (2, 2)-curve 

(x 2 - w 2 )y 2 ± (x 2 - (1 - y 2 )w 2 )z 2 . 

Thus we see that these correspond to the case s = yl — y 2 in the parameterization of 



symmetric (2, 2)-curves given above (4.1), while in Q4.2|) and ( |4.1|) the signs (±) correspond. 

Figures and || display pictures of these four quadrics, together with the x-axis, some 
tangents perpendicular to the x-axis, and the curve on the quadric where the lines are 
tangent. 





(a) (b) 
Figure 1. Real quadrics not meeting the 




(a) (b) 
Figure 2. Real quadrics meeting the x-axis. 

Remark 13. For each of the spheres, there is another sphere of radius r which leads to 
the same envelope, namely the one with center (0, —yo, 0) T . 

The ramification of the (2, 2)-curve of tangents perpendicular to the x-axis is evident 
from Figures H and @. When x = ±1, there is a single tangent line; this line has slope 
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[y, 



[1, 0], i.e., it is a horizontal line. When x = ± 1 — y^ , there is a single tangent 



orem 



12 is evident in these examples. 



line, which is vertical (i.e., which has slope [y, z] = [0, 1]). Figures [I] and |2| depict these 
lines in case they are real. In Figure [I] we have \yo\ > 1, and hence the vertical tangent 
lines are complex. All other values of x give two lines perpendicular to the x-axis and 
tangent to the quadric, but some have imaginary slope. 

The difference in the number of real components of the fiber <yj _1 (C(s)) noted in The- 

The spheres and hyperboloid displayed together 

which interchanges the 



are isomorphic under the change of coordinates z i— > y— T • z, 
transversal tangents of purely imaginary slope for one quadric with the real transversal 
tangents of the other and corresponds to the different signs ± in ( f4.2| ) and ( flip . 

For the sphere of Figure [I], only 4 of the 12 families are real. One consists of ellipsoids, 
including the original sphere, one of hyperboloids of two sheets, and two of hyperboloids of 
one sheet. Since a hyperboloid of two sheets can be seen as an ellipsoid meeting the plane 
at infinity in a conic, we see there are two families of ellipsoids and two of hyperboloids. 
In Figure [j[ we display one quadric from each family (except the family of the sphere), 
together with the original sphere, the x-axis, and the curve on the quadric where the lines 
perpendicular to the tangent to the quadric. 






Figure 3. The other three families. 



Similarly, the hyperboloid of Figure [I] has only 4 of its 12 families real with two families 
of ellipsoids and two of hyperboloids. The sphere of Figure |2] has only 4 of its 12 families 
real, and all 4 contain ellipsoids. In contrast, the hyperboloid of Figure ^ has all 12 of its 
families real, and they contain only hyperboloids of one sheet. 

Many more pictures (in color) are found on the web page^ accompanying this article. 

5. Transversals to two lines and tangents to two spheres 

We solve the original question of configurations of two lines and two spheres for which 
there are infinitely many real transversals to the two lines that are also tangent to both 
spheres. While general quadrics are naturally studied in projective space P 3 , spheres 
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naturally live in (the slightly more restricted) affine space R 3 . As noted in Section 2, 
treat only skew lines. There are two cases to consider. Either the two lines are in R 3 
one lies in the plane at infinity. We work throughout over the real numbers. 



we 
or 



5.1. Lines in affine space R 3 . The complete geometric characterization of configura- 
tions where the lines lie in R 3 is stated in the following theorem and illustrated in Figure |j. 

Theorem 14. Let S\ and S2 be two distinct spheres and let i\ and £ 2 be two skew lines 
in R 3 . There are infinitely many lines that meet i\ and £2 and are tangent to Si and S2 
in exactly the following cases. 

(1) The spheres Si and S2 are tangent to each other at a point p which lies on one line, 
and the second line lies in the common tangent plane to the spheres at the point 
p. The pencil of lines through p that also meet the second line is exactly the set of 
common transversals to £1 and £ 2 that are also tangent to Si and S 2 ■ 

(2) The lines £1 and £2 are each tangent to both Si and S2, and they are images of each 
other under a rotation about the line connecting the centers of Si and S2. If we rotate 
£1 about the line connecting the centers of the spheres, it sweeps out a hyperboloid of 
one sheet. One of its rulings contains £1 and £2, and the lines in the other ruling are 
tangent to Si and S2 and meet £1 and £2, except for those that are parallel to one of 
them. 





(2) 



FIGURE 4. Examples from Theorem 14 



Let £1 and £2 be two skew lines. The class of spheres is not invariant under the set 
of projective linear transformations, but rather under the group generated by rotations, 
translations, and scaling the coordinates. Thus we can assume that 




+ x 



x G 




+ z 




z g 



for some 5 G R \ {0}. As before, there is a one-to-one correspondence between lines 
meeting £1 and £2 and pairs (x, z) G R 2 . The transversal corresponding to a pair (x, z) 
passes through the points (x, 5x, 1) T and (z, —Sz, — 1) T , and has Pliicker coordinates 



[z — x, —5(x + z), —2, —25xz, —(x + z), 5(z — x)) 
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Let S\ have center (a, b, c) and radius r. By Proposition |2| and (O), the transversals 



tangent to Si are parameterized by a curve Ci of degree 4 with equation 

= A5 2 x 2 z 2 + A5(b-a5)x 2 z + ((b-a5) 2 + (l+5 2 )((l+c) 2 - r 2 ))x 2 
(5.1) - 45(b + a5)xz 2 + 2((r 2 -c 2 )(1— 5 2 ) + (1— 6 2 ) + 5 2 (a 2 — l))xz 

-4(l+c)(a+W)x + ((b+a5) 2 + (l+8 2 )((l-c) 2 -r 2 ))z 2 

+ 4(c-l)(a-bS)z + 4(a 2 + b 2 - r 2 ) . 

This is a dehomogenized version of the bihomogeneous equation ( [2.3|) of bidegree (2,2). 
Note also that the curve Ci is defined over our ground field R. The transversals to ii 
and £2 tangent to S2 are parameterized by a similar curve C2. There are infinitely many 
lines which meet £1 and £ 2 and are tangent to Si and S 2 if and only if the curves Ci and 
C 2 have a common component. That is, if and only if the associated polynomials share a 
common factor. We first rule out the case when the curves are irreducible. 

Lemma 15. The curve Ci in ( |5.1| ) determines the sphere Si uniquely. 

Proof. Given the curve ( |5.1|) , we can rescale the equation such that the coefficient of x 2 z 2 
is AS 2 . From the coefficients of x 2 z and xz 2 we can determine a and b, and then from the 
coefficients of x 2 and z 2 we can determine c and r. □ 



Remark 16. By Remark 13, Lemma 15 does not hold if the lines are allowed to live in 



projective space P 3 . We come back to this in Section |5T2 . 



By Lemma [TR there can be infinitely many common transversals to £1 and £2 that are 



tangent to two spheres only if the curves Ci and C 2 are reducible. In particular, this rules 
out cases (1) and (2) of Lemma [| Our classification of factors of (2, 2)-forms in Lemma |3| 
gives the following possibilities for the common irreducible factors (over M.) of Ci and 
C2, up to interchanging x and z. Either the factor is a cubic (the dehomogenization of 
a (2, l)-form), or it is linear in x and z (the dehomogenization of a (1, l)-form), or it is 
linear in x alone (the dehomogenization of a (1, 0)-form). There is the possibility that the 
common factor will be an irreducible (over R) quadratic polynomial in x (coming from a 
(2,0)-form), but then this component will have no real points, and thus contributes no 
common real tangents. 

We rule out the possibility of a common cubic factor, showing that if C\ factors 
as x — xq and a cubic, then the cubic still determines Si. The vector (— 5, — 1, 5xq) t 
is perpendicular to the plane through (xq,5xo,1) t and £ 2 , so the center of Si will be 
(x ,5x , 1) T + A(— 8, — l,5x ) T for some non-zero A £ 1 Thus r 2 = A 2 (l + 5 2 + S 2 xl). 
Substituting this into ( |5.1| ) and dividing by (x — x ) we obtain the equation of the cubic: 

= 5 2 xz 2 + 5{5 2 -l)Xxz + (1+5 2 (1-A 2 ) + 5X(l+5 2 )x )x 

' ,2, +5(\(l+5 2 ) -5x )z 2 + 5(5 2 -l)\x z + 45X + (5 2 \ 2 - 5 2 - l)x . 

Given only this curve, we can rescale its equation so that the coefficient of xz 2 is 5 2 , then 
if 5 7^ ±1, we can uniquely determine A, x and therefore Si, too, from the coefficients of 
xz and x. 
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The uniqueness is still true if 5 = ±1. Assume that 5=1. Then fl5.2Q reduces to 

xz 2 + (2A - x )z 2 + (2 - A 2 + 2Xx )x + 4A + (A 2 - 2)x = . 

Set a := 2A — x , (3 := 2 — A 2 + 2Ax , and 7 := 4A + (A 2 — 2)x . We can solve for A and 
x in terms of a and /5, 

a ± J a 2 + 3/? - 6 -a ± 2 ^a 2 + 3/3 - 6 
= 3 • -0 = § • 

(We take the same sign of the square root in both cases). If we substitute these values 
into the formula for 7, we see that the two possible values of 7 coincide if and only if 
a 2 + 3/5 — 6 = 0, in which case there is only one solution for A and Xq, so a, /5, and 7 
always determine A and xo uniquely and hence Si uniquely. The case 5 = — 1 is similar. 

We now are left only with the cases when C\ and C<i contain a common factor of the 
form x — Xo or xz + sx + tz + u. Suppose the common factor is x — Xq. Then any line 
through p : = {xq, 5xq, 1) T and a point of £2 is tangent to Si. This is only possible if the 
sphere Si is tangent to the plane through p and £2 at the point p. We conclude that if G\ 
and C2 have the common factor x — xo, then the spheres Si and S2 are tangent to each 
other at the point p = (xq, 5xq, 1) t lying on £\ and £2 lies in the common tangent plane 



to the spheres at the point p. This is case (1) of Theorem |14 

Suppose now that C% and C2 have a common irreducible factor xz + sx + tz + u. We can 
solve the equation xz + sx + tz + u = uniquely for z in terms of x for general values of x, 
or for x in terms of z for general values of z, this gives rise to an isomorphism between 
the projectivizations of £\ and £2- The lines connecting q and <p(q) as q runs through the 
points of £\ sweep out a hyperboloid of one sheet. The lines £\ and £2 are contained in one 
ruling, and the lines meeting both of them and tangent to Si are the lines in the other 
ruling. 

Lemma 17. Let H C M 3 be a hyperboloid of one sheet. If all lines in one of its rulings 
are tangent to a sphere S , then H is a hyperboloid of revolution, the center of the sphere 
S is on the axis of rotation and S is tangent to H . 

Proof. We can choose Cartesian co-ordinates such that H has equation x 2 /a 2 + y 2 //3 2 — 
z 2 /^ 2 = 1 for some positive real numbers a, j3, 7. Let the sphere have center (A,B,C) T 
and radius R. The set of points of the form (a, A/5, A7) T , (—a, —A/5, A7) T , (— Xa, /?, A7) 1, 
and (Xa, — /5, A7) 1 " as A runs through K form four lines in one of the rulings. Since the 
two rulings are symmetric, we only need to deal with one of them. 

The sphere S is tangent to a line if and only if the distance of the line from the center 
of S is R. The condition that (A, B, C) T must be at the same distance from the first two 
lines gives the equation 

a(/5 2 + 7 2 )A + ^BC = 0, 

the equality of distances from the othe two lines gives 

/5(a 2 + 7 2 )5 - a^AC = . 

Since a,/5, 7 > 0, the common solutions of these equations have A = B = 0. Using 
this information, the equality of the distances from the first and third lines gives a = j3, 
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or C = ± a/ (a 2 + c 2 )(b 2 + c 2 ) / c. To eliminate this second possibility, consider two more 
lines in the same ruling, the points of the form 

( a T J T ,A7 ) and /?Z ^T' A7 ) 

as A runs through R. The equality of distances from these two lines together with A = 
B = gives a = (3 or C = 0. 

Therefore the only case when (A, B, C) T can be at the same distance from all lines in 
one ruling of H is when a = (3, i.e., if is a hyperboloid of revolution about the z-axis, 
and (A, B, C) T lies on the z-axis. In this case, it is obvious that (A, B, C) T is at the same 
distance from all the lines contained in H, and these lines are tangent to S if and only if 
S is tangent to H. □ 

By this lemma, the hyperboloid swept out by the lines meeting i\ and £ 2 and tangent to 
5*i is a hyperboloid of revolution with the center of Si on the axis of rotation. Furthermore, 
l\ and £ 2 are lines in one the rulings of the hyperboloid, therefore they are images of 
each other under suitable rotation about the axis, the images of i\ sweep out the whole 
hyperboloid, and £1, £ 2 are both tangent to Si. Applying the lemma to S 2 shows that the 
center of 5*2 is also on the axis of rotation and £1, £ 2 are both tangent to S 2 . We cannot 
have 5*i and 5*2 concentric, therefore the axis of rotation is the line through their centers. 



This is exactly case (2) of Theorem [14|, and we have completed its proof. 

5.2. Lines in projective space. We give the complete geometric characterization of 
configurations in real projective space where the line £ 2 lies in the plane at infinity. 

Theorem 18. Let Si and S 2 be two distinct spheres and let £1 lie in M 3 with £ 2 a line at 
infinity skew to £\. There are infinitely many lines that meet £1 and £ 2 and are tangent to 
Si and S 2 in exactly the following cases. 

(1) The spheres Si and S 2 are tangent to each other at a point p which lies on £1, and 
£ 2 is the line at infinity in the common tangent plane to the spheres at the point p. 
The pencil of lines through p that lie in this tangent plane are exactly the common 
transversals to £1 and £ 2 that are also tangent to Si and S 2 . 

(2) Any line meeting £1 and £ 2 is perpendicular to £1 and Si and S 2 are related to each 
other by multiplication by —1 in the directions perpendicular to l\. Thus we are in 
exactly the situation of Remark [73| of Section f| as shown in Figures \^(a) and ^(a). 

Proof. Let IT be any plane passing through a point of £1 and containing £ 2 . Then common 
transversals to £1 and £ 2 are lines meeting £1 that are parallel to II. Choose a Cartesian 
coordinate system in M. 3 such that £1 is the x-axis. Suppose that Si has center (a, b, c) T 
and radius r. Let u = (ui, u 2 , 0) T and v = (vi, 0, t>3) T be vectors with u 2 7^ and v 3 7^ 
parallel to II. Such vectors exist as £1 and £ 2 are skew. A common transversal to £1 
and £ 2 is determined by the intersection point (x, 0, 0) T with £1 and a direction vector 
corresponding to the intersection point with £ 2 , which can be written as u + zv for some 
z G R, unless it is parallel to v. Since Si has at most two tangent lines which meet £1 
that are parallel to v, so by omitting these we are not losing an infinite family of common 
t r ansversals / 1 angent s . 
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The transversals that are tangent to S\ are parametrized by a curve G\ in the xz-plane 
with equation 

= v\x 2 z 2 + u\x 2 + 2v 3 (cvi — av 3 )xz 2 + 2(bu 2 Vi + cuiv 3 )xz 
(5.3) +2u 2 (bu 1 -au 2 )x + ((b 2 + c 2 -r 2 )v 2 -2acv x v 3 + (a 2 + b 2 -r 2 )v\)z 2 

+2((fe 2 + c 2 — r 2 )u\V\ — acuiVs — bu 2 (av\ + cv 3 ))z 
+ ((b 2 + c 2 - r 2 )u\ - 2abu x u 2 + (a 2 + c 2 - r 2 )u\) 

The transversals tangent to S 2 are parametrized by a similar curve C 2 . There are infinitely 
many lines that meet t\ and l 2 and are tangent to S\ and S 2 if and only if C\ and C 2 
have a common non-empty real component. 

It is easy to see from the coefficients of xz 2 , xz and x and the constant term that if 
Ui 7^ or V\ 7^ 0, then C\ determines a, b, c and r 2 and therefore Si uniquely, so if 
C\ is irreducible and «i / or t)i / 0, then there cannot be infinitely many common 
transversals that are tangent to S\ and S 2 . 

Assume now that u± = v\ = 0, this is equivalent to the plane II being perpendicular to 
l\. From the coefficient of x we can determine a, and then from the coefficients of z 2 , z, 
and the constant term we can calculate the quantities a = c 2 — r 2 , (3 = be, and 7 = b 2 — r 2 . 
The equation (a + r 2 )(7 + r 2 ) — (3 2 = is a quadratic equation for r 2 with solutions 

r 2 = - (-a - 7 ± y/(a - 7) 2 + A(3 2 ) . 

Only the larger root is feasible, even when both are positive, since both a + r 2 = c 2 and 
7 + r 2 = b 2 must be non-negative. Hence r 2 , and thus b 2 and c 2 are uniquely determined. 
The values of b 2 , be, and c 2 determine two possible pairs (b, c) which are negatives of 
each other. This is exactly case (2) of the theorem. In fact, this case is illustrated by 
Figures 0(a) and ||(b). 



Let us now consider the cases when C\ is reducible. As in the proof of Theorem [14], we 
need only consider cubics and factors of the form xz + sx + tz + u, x — xq, and z — z . 

Assume that C\ has a component with equation xz + sx + tz + u. As described in the 
proof of Theorem 113, this establishes an isomorphism between the projectivizations of i\ 



and i 2 . The lines connecting the corresponding points of the projectivizations of l\ and 
l 2 sweep out a hyperbolic paraboloid. However, the lines in one ruling of the hyperbolic 
paraboloid cannot all be tangent to a sphere, therefore this case cannot occur. 

Likewise, the factor z — zq cannot appear, since it would mean that all the lines through 
a point of t\ parallel to a certain direction are tangent to S\, which is clearly impossible. 

Consider the case where the equation of C\ has a factor of x — xq. As we saw in the 
proof of Theorem [131 t\ meets the sphere S\ at the point p := (xq,0,0) t , and £ 2 lies in 



the tangent plane to Si at p, and so this tangent plane is parallel to n. 

If x — x is a factor of C 2 , too, then C 2 passes through p and its tangent plane there is 
also parallel to n, so we have case (1) of the theorem. 

To finish the proof we investigate what happens if the common component of C\ and 
C 2 is the cubic obtained from C\ after removing the line x — Xq = 0. 
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The center of Si has coordinates (xo + HU2V3, —fiUiv 3 , —[iU2Vi) T for some \i 6 R, since 
S\ passes through (x , 0, 0) T and its tangent plane there is parallel to II, and we have 
r 2 = ji 2 {u\v1 + u^vf + Substituting this into ( |5.3| ) we obtain the equation of the 

remaining cubic, 

V 3 XZ + U 2 X — V3{XoV3 + 2/iU 2 (f 1 + v 3 ))z 
—AflU^V^Z - U 2 {xqU2 + 2fxv 3 (ul + u\)) = 0. 

If Ui 7^ or v 1 7^ then from the coefficients of this curve we can determine Xq and fi, 
hence ^1 uniquely, so C\ and C2 cannot have a common cubic component. If u% = v± = 
then the above equation factorizes as 

(x - (2fiu 2 v 3 + x ))(vlz 2 + u\) = 0, 

so if C2 contains the curve defined by this equation, then the line x — (2/iw 2 f 3 + xq) = is 
a common component of both C\ and C2, which is a case we have already dealt with. □ 

Appendix A. Calculations from Section || 

We describe the computation of Section |3| in much more detail, giving a commentary 
on the Singular file that accomplishes the computation and displaying its output. The 
input and output are displayed in typewriter font on separate lines and the output begins 
with the Singular comment characters (//). 

The library primdec . lib contains the function sat for saturating ideals, and the option 
redSB forces Singular to work with reduced Grobner (standard) bases. 

LIB "primdec. lib" ; 
option(redSB) ; 

We initialize our ring. 

ring R = 0, (s,t ,a,b,c,d,e,f ,g,h,k,l) , (dp (2), dp(10)); 

The underlying coefficient field has characteristic (so it is Q) and variables s, t, a, . . . ,k, I, 
with a product term order chosen to simplify our analysis of the projection to C 2 , the 
space of parameters. 

We consider the ideal generated by (|3.5| ) 

ideal I = el-g~2, ek-gf, ak-dc, ah-c~2; 

and by the 2x2 minors of the coefficient matrix (|3.6| ). 

matrix M[2] [5] = s , 1-s , -2 , 1-t , t 

al-d"2, 2*(bl-dg), 2* (2bk-cg-df ) , 2*(bh-cf), eh-f~2; 
1=1+ minor(M,2) ; 

We check the dimension and degree (multiplicity) of the variety V(I), first computing a 
Grobner basis for /. 

I = std(I); dim(I), mult(I); 
// 6 8 

Singular gives the dimension of V(I) in affine space C 12 . Since / is homogeneous in the 
variables a,b, . . . ,h,k, I, we consider V(J) to be a subvariety of P 9 x C 2 . Its dimension 
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is one less than that of the corresponding affine variety. Thus V(I) has dimension 5 and 
degree 8. 

In Section ^ we identified three spurious components of V(I) which we remove. The 
first and largest is the ideal of rank 1 quadrics, given by the 2 x 2-minors of the 4x4- 
symmetric matrix ( |3.3p . 

matrix Q [4] [4] = a , b , c , d , 

b , e , f , g , 

c , f , h , k , 

d , g , k , 1 ; 
ideal El = std(minor (Q , 2) ) ; 

We remove this spurious component, computing the quotient ideal (J : Ei). 

I = std (quotient (I, El)) ; dim(I) , mult (I); 
// 5 20 

The other two spurious components describe rank 2 quadrics which are unions of two 
planes with intersection line l\ or £ 2 . 

ideal E2 = g, f , e, d, c, b, a; // intersection line 11 
ideal E3 = 1, k, h, g, f, d, c; // intersection line 12 

The corresponding components are not reduced; rather than take ideal quotients, we 
saturate the ideal / with respect to these spurious ideals. The Singular command sat 
for saturation returns a pair whose first component is a Grobner basis of the saturation 
and the second is the saturation exponent. Here, both saturations have exponent 4. We 
saturate / with respect to E 2 , 

I = sat(I,E2) [1] ; dim(I) , mult (I); 
// 5 10 

and then with respect to E3. 

ideal J = sat (I ,E3) [1] ; dim(J) , mult (J); 
// 4 120 

Thus we now have a variety V(J) of dimension 3 in P 9 x C 2 . We check that it projects 
onto the space C 2 of parameters by eliminating the variables a,b, . . . ,h,k,l from J. 

eliminated, abcdef ghkl) ; 
// _[1]=0 

Since we obtain the zero ideal, the image of V( J) is Zariski dense in C 2 || Chapter 4, §4]. 
However, the projection P 9 x C 2 -» C 2 is a closed map, so the image of V(J) is C 2 . Thus, 
for every smooth (2, 2)-curve C defined by ( |2.5| ), there is a quadric whose transversal 
tangents are described by the curve C. 

We now apply the factorization Grobner basis algorithm facstd to J. The second 
argument of facstd is the list of non-zero constraints which are given by Theorem [?[ 

ideal L = s, t, t-1, s-1, s-t; 
list F = facstd(J,L); 

Singular computes seven factors 
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size (F) ; 
// 7 

Since J and the seven factors Li, . . . , L7 are radical ideals, this factorization can be 
verified by checking that that the following ideals V\ and V 2 coincide. (This part of the 
calculation is archived on the web page''".) 

int i ; 

ideal FF = 1; 

for (i = 1; i <= 7; i++) { FF = intersect (FF,F [i] ) ; } 
ideal VI, V2; 

VI = std(sat(sat(sat(sat(sat(FF,t) [1] ,s) [1] ,t-l) [1] ,s-l) [1] ,s-t) [1]) ; 
V2 = std( sat (sat (sat (sat (sat (J ,t) [1] ,s) [1] ,t-l) [1] ,s-l) [1] ,s-t) [1] ) ; 

Note, in particular, that for any given explicit values of s, t satisfying the nonzero 
conditions, the parametric factorization (in s,t) produced by facstd can be specialized 
to an explicit factorization. 

We examine the ideals in the list F, working over the ring with parameters. 

ring S = (0,s,t), (a,b,c,d,e,f ,g,h,k,l) , lp; short = 0; 

First, the ideal J has dimension 1 and degree 24 over this ring, as claimed. 

ideal JS = std(imap(R, J) ) ; dim(JS) , mult(JS); 
// 2 24 

The first ideal in the list L has dimension 0. 
setring R; FR = F[l]; setring S; 
FS = std(imap(R,FR)) ; dim(FS) , mult(FS); 
// 14 

This ideal is a spurious component from the factorization. It is contained in the spurious 
ideal E 2 . 

FS[5], FS[6], FS[7], FS [8] , FS[9], FS[10], FS[11]; 
// gfedcba 

The other six components each have dimension 1 and degree 4, and each contains a 
homogeneous quadratic polynomial in the variables x and y. 

for (i = 2; i <= 7; i++) { 

setring R; FR = F[i]; setring S; 

FS = std(imap(R,FR)) ; dim(FS) , mult(FS); 

FS[1]; 

print (" ") ; 

} 

// 2 4 

// (-s~2+2*s-l)*k~2+(2*s-2)*k*l+(s*t-l)*l~2 
// 

// 2 4 

// (s-l)*k~2-2*k*l-l~2 



tfrttp : // www . math . umass . edu/~sott ile/pages/212s . html 



COMMON TRANSVERSALS AND TANGENTS IN P 3 



25 



// 

// 2 4 

// (s~2-2*s+l)*k~2+(-2*s+2)*k*l+(-t+l)*l~2 

// 

// 2 4 

// (s~2-2*s+l)*k~2+(-2*s+2)*k*l+(-t+l)*l~2 
// 

// 2 4 

// (s-l)*k-2-2*k*l-l~2 

// 

// 2 4 

// (-s"2+2*s-l)*k"2+(2*s-2)*k*l+(s*t-l)*l"2 
// 

The whole computation takes 7 minutes CPU time on an 800 Mhz Pentium III processor, 
and 3 minutes of that time are spent on the f acstd operation. 

Each of these homogeneous quadratic polynomials factors over Q(\/s, and induces 
a factorization of the corresponding ideal. We describe this factorization — which is carried 
out by hand — in detail for the second component F 2 . We start from the Grobner basis of 
the ideal F 2 computed in the program above, 

(3 - l)k 2 - 2kl - I 2 , (a - l)h + (2t - 2)k + (t - 1)1, fl - gk, 
(A.l) el-g 2 , d + f + g, c, 2b + e, a, 

( s _ l)f k - 2gk - gl, (s - l)f 2 - 2fg - g 2 , ek-fg. 

Over Qiy/s, the first polynomial factors into 

((V^+1)A; + Z) ((yfi-l)k-l) . 

We consider the first factor; the second one can be treated similarly. Substituting I = 
— (-y/5 + l)k into the generator fl — gk, that one factors into 

-k((yfi + l)f + g) . 

Since any zero of F 2 with k = would imply a = c = d = f = g = h = k = l = and thus 
be contained in V(Es), we can divide by k and obtain a linear polynomial. Altogether, 
the first two rows of ( |A.1| ) become a set of seven independent linear polynomials and 
one quadratic polynomial el — g 2 . For any pair (s,t) satisfying ( |2.6[ ) they define a plane 
conic. We leave it to the reader to verify that the three polynomials in the third row are 
contained in the ideal generated by the first two rows. 

In order to show that for none of the parameters s, t satisfying ( p.6|) the map (p vanishes 
identically on this conic, consider the following point p on it: 

(0, -(^i+l)( s -l), 0, -2V£(s-l), 2( v /i+l)( s -l), -2(s-l), 
2(y^+l)(s-l), 4(t-l)-2(t-l)( v /i+l), -2(s-l), 2( v /i+l)(s-l)) 2 . 

The coefficient of w 2 z 2 in <p(C) is 

-4 S (^-l) 2 (v^+l) 2 , 
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so ip(C) does not not vanish identically. 

In order to show that for all parameters s, t satisfying ( |2.6| ) the 12 conies are distinct, 
consider the quadratic polynomials in k and I in the Singular output above. In the 
factorization over Q(y/s, y/i), the ideal of each of the 12 conies contains a generator which 
is linear in k and / and independent of a, . . . , h. To show the distinctness of two conies, 
we distinguish two cases. 

If these linear homogeneous polynomials are distinct (over Q(s,t)), then it can be 
checked that for every given pair (s, t) they define subspaces whose restrictions to (k, I) ^ 
(0, 0) are disjoint. 

In case that the linear homogeneous polynomials coincide then it can be explicitly 
checked that both conies are distinct. For example, both Fi and F§ contain the factor 
(y/s + l)k + I in the first polynomial. As seen above, the corresponding conic of F2 is 
contained in the subspace a = c = 0. Similarly, the corresponding conic of F 5 is contained 
in e = g = 0. Assuming that the two conies are equal for some pair (s,t), the equations 
of the ideals can be used to show further a = b = c= -- - = h = 0. However, due to the 
saturation with the excess component E 2 this is not possible, and hence the two conies 
are distinct. 

The same calculations for the other components are archived at the web page of this 
paper, http: //www. math. umass . edu/~sottile/pages/212s .html. 
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